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This is a study of front dynamics in reaction diffusion systems near Nonequilibrium Ising-Bloch 
ON bifurcations. We find that the relation between front velocity and perturbative factors, such as ex- 

ternal fields and curvature, is typically multivalued. This unusual form allows small perturbations to 
induce dynamic transitions between counter-propagating fronts and nucleate spiral vortices. We use 
these findings to propose explanations for a few numerical and experimental observations including 
O spiral breakup driven by advective fields, and spot splitting. 
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I/") I. INTRODUCTION 

o 

Nonequilibrium Ising-Bloch (NIB) transitions [1] have 
C^Sbeen identified recently as important mechanisms of 
^^pattern formation in reaction-diffusion [2,3] and liquid 
^^crystal [4] systems. Mathematically, a NIB transition 
'"T^amounts to a pitchfork bifurcation where a station- 
^ary (Ising) front loses stability to a pair of counter- 
1 1 .p ropagating (Bloch) fronts. The coexistence of two Bloch 
"^fronts far beyond the bifurcation allows the formation of 
(^r egular patterns such as periodic traveling domains and 
rotating spiral waves. 

In this paper we show that near a NIB bifurcation 
spontaneous transitions between the two Bloch fronts be- 
come feasible and that a number of apparently unrelated 
phenomena can be understood in terms of these transi- 
tions. Front transitions of this kind can be induced by 
extrinsic perturbations, such as advective fields, or in- 
trinsic perturbations, such as front interactions and cur- 
vature. They may occur uniformly along the front, re- 
versing its direction of propagation, or locally, nucleat- 
ing spiral- vortex pairs. We study the mechanism of such 
transitions and suggest explanations for two experimen- 
tal observations: spiral breakup induced by the onset of 
convection in chemical reactions [5], and spot splitting 
[6]. We suggest that these experiments, the transition 
from labyrinthine patterns to spiral turbulence found nu- 
merically in [7], and possibly other breakup phenomena 
[8,9] are all realizations of the same mechanism. 
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Perturbation 

FIG. 1. A typical front- velocity vs. perturbation graph 
near a NIB bifurcation. The figure is a plot of V vs. I ac- 
cording to the stationary form of (13) with cto = — 2.0 and 
ai = 5.0. 



The key feature underlying front transitions is a mul- 
tivalued dependence of front velocity on extrinsic and 
intrinsic perturbations close to a NIB bifurcation. Fig. 
1 shows a typical form of such a relation. The up- 
per and lower branches correspond to the two counter- 
propagating Bloch fronts. At least one of them termi- 
nates at a realizable perturbation strength. A perturba- 
tion that drives the system to the end of a given Bloch 
front branch induces a transition to the other branch as 
indicated by the arrow. 
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We consider an activator-inhibitor reaction-diffusion 
system exhibiting a NIB bifurcation [3]: 



vV 



v t = e(u — a\v — do) + J • Vi> + 6V 2 v, 



(la) 



(lb) 



where J represents an advective vector field. In the con- 
text of chemical reactions u and v represent concentra- 
tions of two key chemical species, and J may stand for 
a convective flow field. We note that there is no a pri- 
ori reason to assume that both species are advected dif- 
ferently by the flow field. However, chemical reactions 
normally involve many more than two chemical species. 
Two-variable models like (1) are obtained by adiabatic 
elimination of fast reacting species, and in this process 
transport terms are renormalized. As shown by Dawson 
et al. [10] this can lead to an effective differential flow as 
assumed in (1). The external field J may also stand for 
an electric field which advects ionic species [8,11]. 

In section 2 we use equations (1), simplified by set- 
ting 8 = 0, to derive an evolution equation for the front 
velocity near a NIB bifurcation. We obtain a multival- 
ued relation between the front velocity and the advective 
perturbation J, and show that transitions between the 
different front velocity branches follow from the gradient 
nature of the evolution equation. In section 3 we study 
the relation between the front velocity and the advec- 
tive perturbation for 8/e 1 and at arbitrary distance 
from the NIB bifurcation. The multivalued nature of 
this relation near the NIB bifurcation is used to demon- 
strate how spiral breakup can be induced by an advective 
field. The relation to convection induced disorder in the 
Belousov-Zhabotinsky (BZ) reaction [5] is discussed. In 
section 4 we show that multivalued relations between the 
front velocity and intrinsic curvature perturbations near 
a NIB bifurcation can be used to explain spot splitting [6] 
and transitions from labyrinthine patterns to spiral tur- 
bulence [7]. We conclude in section 5 with a discussion 
of further possible implications of this work. 



a summary of the main findings. The system (1) has a 
stationary Ising front solution, 



u (x) 



tanh(t/_|_a;/V2), vq(x) = a 1 1 uq(x), (2) 



connecting the up state at x = — oo to the down state 
at x = oo. This front solution exists for all e > 0. For 
e < e c = aj~ 2 two additional Bloch front solutions appear, 
propagating at velocities 



±[^T(e c -e)]^. 



(3) 



They correspond to an up state invading the down state 
(positive speed) and to a down state invading the up state 
(negative speed). Their leading-order forms for |c| <C 1 
are 

u(x, t) = uo(x — ct), v(x,t) = a^ 1 uo(x — ct + ca\), (A) 

[3]. To that order, the Bloch front structure differs from 
the Ising front structure in that the v field is translated 
with respect to the u field by an amount proportional to 
c. 

Our objective is to derive an evolution equation for 
the front velocity, C(t), near the NIB bifurcation, whose 
solutions describe dynamic transitions between the two 
Bloch fronts. We introduce a normalized front velocity, 
V(t), such that V = corresponds to the Ising front and 
V = ±1 to the Bloch fronts. The actual front velocity 
is given by C = cV , where c is a positive constant to 
be identified with the absolute value of the Bloch front 
velocity given by (3). We derive an evolution equation 
for V using asymptotic expansions in the small positive 
parameter c and demanding uniform validity. We write 

oo 

u(z,T) = u (z) + y^c n u n (z,T), (5a) 



■j(z,T) = vo(z) + J2 c "M^T), (56) 



II. THE DYNAMICS OF FRONT TRANSITIONS 



(5c) 



To study front transitions we consider the simpler case 
of a non-diffusive v field, 6 = 0, in one space dimen- 
sion, x, and assume first a symmetric system, ao = 0, 
with no external field, J = 0. We choose a\ such 
that equations (1) describe a bistable medium having 
two linearly stable, stationary homogeneous states: an 
up state, (u + ,v + ), and a down state, («_,«_), where 



Va-ar 1 ). 



and v 



+ 



In Ref. [3] we studied one-dimensional front solutions 
of (1) propagating at constant speeds. The following is 



where z = x — C(T) and T = c 2 t. In terms of the front 
position, C, the normalized velocity is given by V(T) = 
c _1 C t = c(t ■ We use these expansions in (1) and solve 
for the corrections, u\, «2, ... and v\, i>2, to the fields. 
At order c we obtain 



where 



- 1 - 



-V 



3«o 



UOz 
VOz 



-1 

-e ai 



(6) 



(7) 
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Since the adjoint operator, C\ has the null vector 
(uo z , —e ~ 1 vo z ) T , a solvability condition leads to eo = 
e c = aj~ 2 . Solving (6) we find 



«i = 0, v 1 = Vu 0z - 
At order c 2 we obtain 

«2 \ ( 



C 



V 2 



-eiai«i — V^i 



(8) 



(9) 



Solvability of (9) yields e\ = and the solutions 

u 2 = ^a 1 V 2 zu 0z , v 2 = ^V 2 zu 0z + a 1 V 2 u 0zz . (10) 



Finally, at order c 3 



«it - Vv 2z + aie 2 vi 



(11) 



Solvability of (11) gives the evolution equation for the 
front velocity; 



V T = vV(l - V 2 ), 



e < e. 



(12) 



where fi = ai\e 2 \ and e 2 = u 0z u 0zzz dz / jT^Kz = 
— 2m^_/5. Equation (12) reproduces the Ising and the 
Bloch front solutions of (1) as constant speed solutions 
with V = and V = ±1, respectively. Since e = e c — 
|f2|c 2 +0(c 4 ), the Bloch front velocities, C = ±c, coincide 
with (3). Equation (12) also contains information about 
the stability of these solutions (for e < e c ) with respect 
to Galilean boosts; the Ising front is unstable and the 
Bloch fronts are stable. 

Next, we consider the nonsymmetric case, do ^ 0, with 
constant advective perturbation J = Jx. For simplicity 
we take do = c 3 ao and J = c 3 7, where c*o and I are of 
order unity. The order parameter equation now reads 

Vt = ftV(l — V 2 ) + veto + —I, (13) 

ai 



where v 



One immediate result from (13) 



\/2 a i ( a i — 1 ) 

is a multivalued relation, of the form shown in Fig. 1, 
between the constant velocity of a front and the advec- 
tive perturbation I (or J). Additionally, equation (13) 
describes a gradient flow (despite the nongradient nature 
of (1)) derivable from the potential 



n V2 \ V2 



1 



veto ■ 



In a range of I values, I m i n (a ) < I < I max (a ), U{V) 
is a double well potential. The wells at V = V_ < 
and V = V+ > pertain to Bloch fronts propagating 
to the left and to the right, respectively. Decreasing I 
below I m in causes the well at V+ , or the Bloch front that 
propagates to the right, to vanish. Since the dynamics are 



gradient, convergence to the well at V = V_ follows. This 
amounts to a transition from a Bloch front propagating 
to the right (upper branch in Fig. 1) to a Bloch front 
propagating to the left (lower branch in Fig. 1), induced 
by the advective perturbation J. 

We have considered the case 6 = 0, which is simpler for 
analysis, and showed that a perturbation driving a given 
Bloch front branch past its end point induces a transi- 
tion to the other Bloch front branch. In the following 
sections we consider the more general case, 8 ^ 0, by 
studying constant speed front solutions for parameters 
satisfying 6/e 1. The analysis culminates in velocity - 
perturbation relations valid both near and far from the 
NIB bifurcation, but does not contain the information 
about the dynamics of front transitions derived in this 
section. Numerical results, however, indicate that the 
same dynamical behavior holds for 6 ^ as well; a front 
transition is induced when a given front reaches the end 
point of a velocity branch. 



III. SPIRAL BREAKUP INDUCED BY AN 
ADVECTIVE FIELD 

The earliest observations of spiral breakup in chem- 
ical reactions occurred in open cells containing the BZ 
reagents [5]. Convection in the form of Benard cells, 
induced by evaporative cooling of the free surface, ini- 
tiated processes of spiral breakup and spiral wave nucle- 
ation that destroyed the order of the pre-existing chemi- 
cal pattern. We do not attempt a complete description of 
this system. Instead we use a simplified model, including 
what we believe are the essential ingredients responsible 
for this behavior (see discussion at the end of the section). 
We consider equations (1) in two space dimensions near 
the NIB bifurcation, assume 6/e 1, a condition nor- 
mally met in experiments on the BZ reaction, and take a 
time independent advective field with hexagonal spatial 
structure (to model the convective Benard cells). 

Before simulating equations (1) with an hexagonal ad- 
vective field, we derive the relation between the front 
velocity and a constant advective field J for 6/e ^> 1. 
This relation will help us to identify conditions leading 
to spiral breakup. In Refs. [3,12] we studied the NIB bi- 
furcation for 6/e ^> 1 in the absence of an advective field 
(J = 0), using a singular perturbation approach (see also 
Ref. [13]). In this parameter range the spatial variation 
of the v field is on a scale much longer than that of u. In 
the narrow front region where the variation of u is of or- 
der unity, v is approximately constant. Solving (1) (with 
J = 0) in this inner region yields the relation 



Co = 7=v f , 



_3_ 



between the front velocity, Co, and the (small) constant 
value, Vf, of the v field at the front position. Away from 
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the narrow front region the u field is enslaved to the v 



field, u 



u±(v), where u±(v) solve u 



0. 



Solving (1) (with J = 0) in the outer regions on both 
sides of the front and matching the solutions at the front 
position yield another relation between Co and Vf, 



Co 



a 



q 2 (C 2 + 4r] 2 q 2 y/ 2 q 2 



Eliminating Vf we get the following implicit relation for 
the front velocity [3,12]: 



C = F(C ,r ] ), 



F(X,Y) 



3X 



V / 2q 2 (X 2 +4q 2 Y 2 y/ 2 
eh, q 2 = ai + 1/2, and C a 



3a n 
y/2q 2 



(14) 



A 



where r/ 2 

graph of Co vs. rj (or Co vs. e at constant 6) yields a 
NIB bifurcation diagram valid at any distance from the 
bifurcation point as long as 6/e 1. 

Equations (14) can be used to find a relation for the 
front velocity C in the presence of a constant pertur- 
bation J. Consider a planar front solution of (1), with 
J = Jx, propagating at a constant velocity C in the x 
direction. It satisfies the equations 



Ci 



0, 



6v" + (C + J)v' + e(u - aiv - a ) = 0, 



(15a) 



(156) 



6jv" + Cv' + €j(u - aiv - a ) = 0, 



where the prime denotes differentiation with respect to 
the single independent variable x—Ct. Let us now rewrite 
equations (15) in the form 

u" + Cu' + u-u 3 -v = 0, (16a) 

(166) 

\j, oj = 8Aj, and Aj = ^y. According 
to (16), front solutions of (15), the perturbed system with 
parameters e and 6, are equivalent to front solutions of 
an unperturbed system with effective parameters ej and 
6j. Using relation (14), obtained for the unperturbed 
system, and replacing Co, e, and 6 by C, ej, and 6j, 
respectively, we find 



where ej = eAj, 6j 



C = F(C + ,J, ry), 



(17) 



where F is given by (14). Figs. 2a — 2c show graphs of C 
vs. J, obtained by solving (17) for C deep in the Bloch 
regime (2a), near the NIB bifurcation (26), and deep in 
the Ising regime (2c). The single velocity branch deep in 
the Ising regime folds into three connected branches near 
the NIB bifurcation with a multivalued form resembling 
the V vs. I graph obtained from (13) for 6 = 0. Deep in 
the Bloch regime the three branches become practically 
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FIG. 2. Front velocity, C, vs. advective perturbation, J, 
and curvature, K, for 8/e 3> 1: (a,d) deep in the Bloch 
regime, (b, e) near the NIB bifurcation, (c, /) deep in the Ising 
regime. 

disconnected; small J perturbations cannot induce front 
transitions. 

We now demonstrate how an hexagonal advective field, 
modeling the Benard convective cells, can induce local 
front transitions leading to spiral breakup. We numeri- 
cally integrated (1) with Neumann boundary conditions 
and J = V</> where <f> = A^2 i=1 cos(q 8 - • x) with qi = 
(Q,0), q 2 = (-Q/2.V5Q/2), q 3 = (-0/2,-^30/2). 
We chose parameters e and 6 to obtain a C — J rela- 
tion as depicted in Fig. 3. The upper velocity branch 
terminates at J = J+ < 0, and the lower branch at 
J = J~ > 0. The peak value of the advective field, 
J max = Max{|J|}, was chosen such that 

I c I ^ Jmax ^ Jc • 

With this choice of parameters the advective field may 
induce front transitions from the upper velocity branch to 
the lower one but not vice versa. The initial conditions, 
shown in Fig. 4a, pertain to a regular spiral wave in the 
absence of the advective perturbation. The dark area 
represents an up state domain. The leading front of the 
spiral wave (an up state domain invading a down state) 
and the trailing front (a down state invading an up state) 
correspond to the upper and lower velocity branches in 
Fig. 3, respectively. 

Figs. 46, c, d show the time evolution of the initial spi- 
ral pattern. As the spiral arm propagates through the 
hexagonal advective pattern, alternating portions of its 
leading front experience phases during which J = J • n < 
0, i.e. the advective field has a component pointing in 



4 



= 1 




-0.3 j; j- 0.3 

J 

FIG. 3. The relation between the front velocity and the 
advective field for the parameters used in the spiral breakup 
simulation shown in Figs. 4 and 5. 

a direction opposite to the direction of propagation, n. 
With the above choices of parameter values the advec- 
tive field is strong enough to drive portions of the leading 
front from the upper branch, past its end point, to the 
lower branch. These local transitions nucleate a spiral- 
vortex pair for each hexagon (Benard cell) the stripe 
crosses. These nucleation processes are evident in Figs. 
5 (crossings of the u and v zero contour lines) which pro- 
vide closer looks at the processes occurring in Figs. 4. 
Portions of the leading front that undergo front transi- 
tions eventually annihilate with the trailing front, break- 
ing the stripe into disjoint pieces. Note that the trailing 
front does not undergo front transitions because the ad- 
vective field is not strong enough to drive fronts on the 
lower velocity branch to its end point. 

The actual experimental situation is more complex. 
The system is three dimensional (the flow at the bot- 
tom is in the opposite direction to the flow at the top) 
and we have ignored possible feedback effects of the re- 
action on the advective field. The mechanism described 
above can therefore account only for the initial spiral 
breakup behavior and not for the asymptotic state. We 
have also considered a bistable medium while the exper- 
imental system is excitable, but we expect excitable sys- 
tems to exhibit similar qualitative behaviors. The front 
structures that appear in excitable systems coincide with 
those in bistable systems, to leading order in e <C 1, and 
we may expect them to undergo an instability analogous 
to the NIB bifurcation. Since single front structures do 
not exist in excitable systems we should rather look at 
the behavior of front pairs or pulses. Indeed, the NIB bi- 
furcation in bistable systems leads to a pulse instability 
that has a close analog in excitable systems; depending 
on the value of 6, traveling pulses in both systems either 
collapse or develop into breathers as e is increased past 
a critical value [3,14,15]. 




FIG. 4. Breakup of a spiral wave induced by a hexagonal 
advective pattern. The light and dark regions correspond to 
down and up states, respectively. The dotted curves denote 
contours of constant advection speed. The convection flow 
direction is outward from the centers of the hexagons. Frame 
a is the unperturbed spiral wave and frames b, c, d are taken 
at times t = 100, 140, 220 from the onset of the advective 
pattern. Parameters used: cio = —0.1, cii = 2.0, e = 0.032, 
8 = 0.9, A = 1.59 and Q = 0.06283. 

IV. VORTEX NUCLEATION AND DOMAIN 
SPLITTING INDUCED BY CURVATURE 

So far we studied the effects of an extrinsic perturba- 
tion in the form of an advective field. Even without ex- 
trinsic perturbations there exist intrinsic factors, such as 
curvature and front interactions, that affect planar front 
propagation. We now address the effect of curvature near 
a NIB bifurcation. In Ref. [7] we have found the follow- 
ing relation between the normal velocity of a front, C , 
and its curvature, K: 

C = F(C + 6K, rj) - K, (18) 

where F is given by (14). Graphs of C vs. K deep 
in the Bloch regime, near a NIB bifurcation, and deep 
in the Ising regime, are shown in Figs. 2d, 2e and 2/, 
respectively. A comparison of these figures with Figs. 2a- 
2c indicates that the multivalued dependence depicted in 
Fig. 1 is a general feature of front dynamics near a NIB 
bifurcation. Perturbations of different nature, curvature 
and an advective field in this case, can have the same 
effect of inducing front transitions. 

Recently, an intriguing pattern formation phenomenon 
has been observed in numerical and laboratory experi- 
ments [6,16]. An expanding chemical spot was found to 
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FIG. 5. A closer look at a typical breakup process in Fig. 
4. The thick (thin) lines are contours of u = (v = 0). 
The direction of front propagation follows from the rule that 
the v = contour always lags behind the u = contour. The 
frames a, b, c, d pertain to times t = 140, 160, 180, 200. They 
show a local front transition, accompanied by the nucleation 
of a vortex pair (the crossing points of the contour lines), and 
the breakup of the up state domain. Parameters are the same 
as in Fig. 4. 



lose its circular shape and split into two spots. New- 
born spots followed the same course of events eventually 
filling the reaction domain. We suggest that spot or do- 
main splitting can result from dynamic front transitions 
induced by curvature variations. Imagine a C vs. K re- 
lation in which the upper branch terminates at a small 
positive curvature value, A'+ > 0, as depicted in Fig. 6. 
The leading front of an expanding disk-like domain cor- 
responds to the upper branch (an up state invading a 
down state). As the domain expands the front curvature 
decreases. When it falls below, A'+, the critical curva- 
ture pertaining to the edge of the upper branch, a front 
transition occurs. As a result the domain stops expand- 
ing and begins shrinking. We have verified this scenario 
numerically using a circularly symmetric version of (1) 
(with J = 0). Two-dimensional realizations of disk-like 
domains are never perfect; there are always flatter parts 
of the front which undergo the transition first. Like in 
the case of advective perturbation, these local front tran- 
sitions nucleate spiral vortex pairs and lead to domain 
splitting as shown in Fig. 7 for an initially oval-shaped 
domain. 

The transition from labyrinthine patterns to spiral 
turbulence found in [7] follows from a similar velocity- 
curvature relation except that the upper branch termi- 
nates near zero or at negative curvature values. In that 



c 




-0.3 



0.3 



FIG. 6. The relation between front velocity and curvature 
for the parameters used in the spot splitting simulation shown 
in Figs. 7. 




a 







FIG. 7. Spot splitting induced by curvature variations. 
The frames a, b, c, d pertain to times t = 80, 240, 280, 340. 
Local front transitions occur at the natter portions of the 
expanding front. They are accompanied by nucleation of vor- 
tex pairs, and followed by spot splitting. Parameters used: 
ao = -0.15, ai = 2.0, e = 0.014, 8 = 3.5. 
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case the mechanism that drives the leading front past 
the end point is the transverse instability which produces 
segments with negative curvature. 

V. CONCLUSION 

We have shown a few examples of local front transi- 
tions induced by extrinsic or intrinsic perturbations near 
a NIB bifurcation. Such transitions are accompanied by 
the nucleation of spiral vortex pairs, bounding the front 
segments that underwent the transitions. The spiral vor- 
tices can be viewed as "front" structures along the front 
line where the v field goes from a positive value (pertain- 
ing to one Bloch front) to a negative value (pertaining to 
the other Bloch front), and vice versa. The nucleation of 
spiral vortices is followed by domain splitting or breakup 
as Figs. 5c, d and 76, c, d show. We note that domain 
breakup is not a necessary outcome of spiral vortex nu- 
cleation. It occurs when the speed at which two fronts 
approach one another is high enough relative to the rate 
of diffusive dissipation of v in the region bounded by the 
fronts. Then, the fronts can approach one another to 
within a distance of order unity where diffusive dissipa- 
tion of u becomes effective and can annihilate the fronts. 
This is usually the case near the NIB bifurcation and 
beyond it (i.e. in the Bloch regime). 

The ideas presented in this paper may explain many 
more observations. Most recently Taboada et al. re- 
ported experimental observations, in the BZ reaction, of 
circular wave breakup phenomena induced by an electric 
field [8]. They also studied the two- variable Oregona- 
tor model with an advective term and numerically repro- 
duced the experimental results. Among their findings is 
a monotonically increasing relation between the critical 
electric current amplitude (the advective field) needed 
for wave breakup and the "excitability", 1/e. This find- 
ing can be understood once we consider the change in 
the C — J relation near the NIB bifurcation as 1/e is 
increased. As Figs. 2a,b,c imply, the critical J value, 
J^"(e) < 0, at which the upper velocity branch termi- 
nates, increases in absolute value as we go farther into the 
Bloch regime or increase 1/e. In other words, a stronger 
field J is needed to cause front transitions (and conse- 
quently breakup) as the excitability (1/e) is increased. 

The onset of breathing motion in pulses near the NIB 
bifurcation [3,14,17] can also possibly be interpreted in 
terms of dynamic front transitions. In this case the tran- 
sitions are induced by intrinsic front interactions that 
become significant as two fronts approach one another. 
A multivalued front velocity relation as depicted in Fig. 
6 with the horizontal axis replaced by the reciprocal dis- 
tance between the fronts can account for a breathing mo- 
tion. Such a relation has not been derived yet. 

We have studied an activator-inhibitor type reaction 
diffusion model but we expect the basic results to be ap- 



plicable to other systems undergoing a NIB bifurcation. 
These include periodically forced oscillatory reactions [1] 
and liquid crystal systems [4]. 
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